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' For the S-states of positronium and muonium, the terms of an expansion of 

energy levels in powers of the fine structure constant a are also members of a 
"recoil series". The first two terms of that series are calculated to all orders 
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The calculation of energy levels to the order a 6 for the S-states of positronium and 
muonium (e~/i + ) by Pachucki JTJ-Q has recently been confirmed for equal masses by en- 
tirely analytic methods Both calculations use non-relativistic quantum electrodynamics 
(NRQED) ||. The only two-fermion equation which is solved nonperturbatively is the 
Schrodinger equation for a particle of reduced mass /x = 11111112/ 'm (m = ni\ + 1712) in a 
Coulomb potential V = —a/r (H = c = 1). Apart from a state-independent function 
F{fi/m) which is only numerically known for arbitrary 1711/1712, all terms to order a 6 and 
a 6 log a can also be arranged in a finite series in fi/m for the total cms energy E, 

E = m + E l + E 2 + E log + E 3 , (1) 

where E{ is of the order /i(/i/m) l_1 , and E\ og contains logarithms of mass ratios. In the 
following, Ei and E 2 are derived to all orders in a, and most of the sixth-order terms 
E^ ~ a 6 /i 3 /m 2 of E^ are confirmed (radiative corrections will be omitted). In combination 
with the complete expression for E^\ the series (P increases the precision of QED bound 
state calculations for systems with mi 7^ vri2- Our method also reduces the gap between 
relativistic and non-relativistic expansions. 



The history of relativistic recoil expansions is long and sad Frequently, one starts 

from a Dirac equation for an electron of mass mi in the potential V(r). The subsequent 
evaluation of recoil corrections of the order of mi/m 2 shows that most terms can be taken 
care of by replacing mi by [i in the Dirac equation. Empirically then || , Ei follows from that 
equation. Similarly, the /i 2 /m-form has been found to order a 4 for the hyperfme splitting, 
which is part of E 2 . But there has been no indication that to a given order of a, the series 
(|I|) would end after a finite number of terms. For some P-states, hyperfme mixing requires in 
fact infinitely many terms already at the order a 4 . The quantitative mixing disagrees with 
the standard hyperfme operator of the Dirac equation. On the contrary, this mixing follows 
easily from the Breit operators of NRQED. At the end, the Dirac equation approach has been 
completely abandoned for two-body systems. It is then remarkable that the nonperturbative 
use of the Schrodinger equation leads to the form (Q), in which Ei is given precisely by the 
Dirac equation with reduced mass. For a state of angular momentum j, principal quantum 
number n and with the abbreviation j + ~ = j + , the a 6 -part of E\ is 

Ef ] = -/i« 6 (V + 3n~% 2 - 6n-%> + V 3 ) /8n 3 , (2) 
from which Pachucki's result follows for j + = 1. 

We have recently derived a relativistic, Dirac-like two-fermion equation from perturbative 
QED ||, which explains this mystery and renders the calculation of E 2 trivial, again to 
all orders in a. The evaluation of E\ og and of the state-independent F(p/m) requires the 
calculation of loop integrals, which remains to be done. The progress arises from the strict 
use of relativistic two-body kinematics in the first Born approximation, and from the reduc- 
tion to 8 x 8 components before Fourier transforming. Any system of two free particles of 
masses m x and m 2 satisfies the cms equation (p 2 — k 2 )ip = 0, where the eigenvalue k 2 can 
be expressed in terms of E 2 , m\ and m\. It can be brought into the form k 2 = e 2 — 

p E = mim 2 /E, e=(E 2 -m\- m 2 2 ) /2E. (3) 

The free equation is converted into an explicit eigenvalue equation for E 2 by the substitution 

p = fi E r, p p = p/n E , e = e/fi E = {E 2 - m\ - ml)/2m 1 m 2 . (4) 

When the interaction is added, the Coulomb potential V(r) is transformed into V(p). The 
resulting dimensionless Dirac equation is 

((3 + l5 ((Ti + o- hf )p p + V(p)-e)^ = 0, (5) 



cr hf = -i(Ti x er 2 V(p)m 1 m 2 /E 2 . (6) 

The &1 and a 2 are Pauli matrices; the product 75^! is normally written as ex. For compar- 
ison, the hyperfme cr h f of atomic theory sets E = m 2 and replaces Vp by [V,p]/2. 
To begin with, we evaluate the hyperfme energies E^f by first-order perturbation theory. 
For orbital angular momentum I — j ± | and total angular momentum / = j ± |, they are 

E hf = 2^a 42 4^4(jle-K D /2)/(^ - 7), (7) 
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K D = 2{l-j)j + , 1 2 =f + -a 2 . (8) 

Special cases of this formula are found in |T0]. To order a 6 and for k d = the quotient 
of the last two brackets in (|?p is 



] 2 + e-K D /2_ 1 



47 - 7 2jj_ 



1 + a 2 
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+ 2f + 2n 3+ 2n 2 2n 2 (j + 1) 



(9) 



The hyperfine splitting is defined as AE h f = E(f — j + ~) — E(f = j — |). For S-states, 
I = implies j + = 1, 2(f — + = 1 ^ or / = anc ^ — ^ ^ or / = ^> which makes a factor 
I in AEhf. Pachucki's result for the part E^lf °f follows by approximating E 2 pa m 2 
in the denominator of (0). The remaining two terms of E 2 appear in the hyperfine-averaged 
shift E = \E(f = 1) + \E(f = 0). They are conveniently evaluated by a non-relativistic 
reduction: 

With the approximation E 2 pa m 2 in the hyperfine operator, (H) is an explicit eigenvalue 
equation that permits the standard reduction by elimination of the small components. The 
resulting Schrodinger equation is 

(1 + p 2 p /2 + V(p) - esMsch = 0. (10) 

It has the familiar non-relativistic eigenvalues, e,s c h — 1 = —a 2 /2n 2 . The equation becomes 
quite powerful when its centrifugal barrier 1(1 + l)/p 2 is replaced by an effective barrier 
l'(V + l)/p 2 , which includes the lowest-order spin-orbit and hyperfine couplings: 

V 2j + f + l m) 1 } 

The principal quantum number n = n r + I + 1 gets replaced by n* = n r + 1 1 + 1 = n + 51, 
and the eigenvalues Eg ch follow from (|l(J) as 

Eg ch — m 2 = —a 2 mim2/n* 2 pa -« 2 mim2(l — 5l/n) 2 /n 2 . (12) 

To order a 6 , one obtains 

E Sch -m= -a 2 (l - 5l/n) 2 fi/2n 2 - «V(1 - 45l/n)/8mn 4 - a 6 /i 3 /16m 2 n 6 . (13) 

This expression contains both terms of E^ and two of the three terms in E^ |2| (note that 
(SI) 2 is quadratic in the hyperfine interaction). The third term arises from the S-D-mixing 
in second order perturbation theory and is not calculated here, Egjj = — |5/ i ia 6 /i 3 /m 2 n 5 . It 
appears is E 3 with a factor |, and in AE 3) hf with a factor 1. 

There are two more n _5 -contributions to AE^f, one from the hyperfine part of 51, Shf = 
2a 2 -^(f —j) /(/ + §) in fll~3l ) (in the bracket following a A n 2 ), and one from setting E 2 = m 2 — 
a 2 rriim2/n 2 in (|7|), which effectively enlarges all hyperfine effects by a factor 1 + a 2 fj,/mn 2 . 
The total coefficient of a 6 /i 3 /m 2 n 5 is then — | + | + | = 3 j 2 -, in agreement with [p]]. 
For higher orders in a, the expansion ([I]) is conveniently replaced by the simpler expansion 
for e — 1 = (E 2 — m 2 )/2m 1 m 2 

e = 1 + ix + e 2 + S\ og + £3, (14) 
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in which e\ is pure "Dirac" as in E±, and e 2 is pure "hyperfine". In other words, the non- 
hyperfine terms of E 2 are canceled to all orders in a. 

There is one more hyperfine contribution to AE^hf [|I| which does not follow from (0): 

^^^fhl + i-J-^). (15) 

It indicates that the combination of Dirac and hyperfine operators is still incomplete, at least 
for S-states. This is seen also for mi = m%, where the combination (3 + V — e of (|5|) reduces 
to V — E 2 /2m\ in the small components. In parapositronium, the same factor appears in the 
combination (cri + (Thf)p in the large components and cancels out. The resulting differential 
equation is reduced to the confluent hypergeometric one, in a shifted variable r' = r — 2a/ E. 
Its singularity at r' — occurs at r = 2a/ E, which excludes any bound state interpretation. 
Further progress in the calculation of E3 and E\ og in ([!]) will require the recalculation of 
Feynman diagrams without kinematical approximations. The resulting formulas should not 



depend on the sign of E, as a consequence of CP-invariance In the integration over 

loop momenta, one normally expands the zero- components p® of the fermion momenta about 
their individual external values. To maintain CP-invariance the product S1E2 of the fermion 
spin summations = pcii + m« may be rearranged as follows: 

EiS 2 = 2pi7ip 2 72 + 2mim 2 - (ptji - m 1 )(p 2 j2 - m 2 ). (16) 

The last product contributes little; it vanishes on either mass shell, but its CP-invariant 
evaluation would have to include negative values of p\ and p®- In the first term, on the other 
hand, one may keep p\p\ positive without breaking CV. 

The appearance of a factor E~ 2 in the hyperfine operator (P) has been confirmed here for the 
first time, by comparison with the a 6 /i 3 /m 2 -terms [0]. Outside QED, a similar effect seems 
to exist in the hyperfine splitting between vector (1~) and pseudoscalar (0~) mesons. The 
splitting in E 2 , A = E 2 (l~) — E 2 (0~), increases uniformly from 0.48 GeV 2 for the heavy 



b quarkonium |L2| to 0.57 GeV 2 ~ m 2 for the p — ix splitting [|13|]. If the light pseudoscalar 
mesons could ever be treated by a potential model, their E 2 (0~) would be lowered by the 
occurrence of 1/E 2 in their hyperfine operator. The effect would become extreme in the 
limit E 2 — > 0. 
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